Distribution of Non-uniform Demagnetization Fields in Paramagnetic Bulk Solids 
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A general calculation for the distribution of non-uniform demagnetization fields in paramagnetic 
bulk solids is described and the fields for various sample geometries are calculated. Cones, ellipsoids, 
paraboloids and hyperboloids with similar sample aspect ratios are considered. Significant differ- 
ences in their demagnetization fields are observed. The calculation shows that the demagnetization 
field magnitudes decrease along the axis of symmetry (along z) where an externally applied magnetic 
field is aligned, and increase in the vicinity of the lateral surfaces with the largest field values found 
in the cone and the narrowest field distributions found in the hyperboloid. Application is made to 
the theoretical modeling of the ^H-NMR spectra of a single crystal of field-induced superconductor 
A-(BETS)2FeCl4 with a rectangular sample geometry, providing a good fit to the measured NMR 
spectra. This calculation is also applicable to diamagnetic or ferromagnetic materials in general. 

PACS numbers: 75.30.Cr, 75.20.-g, 75.30.-m, 76.60.-k 
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I. INTRODUCTION 

The demagnetization field as a type of local magnetic 
field contributes to both an inhomogeneous broadening of 
nuclear magnetic resonance (NMR) spectra and a shift in 
the NMR resonance frequencies. Similar effects could be 
produced by other local field sources at a nucleus as well 
when they are not negligible. Thus demagnetization 
fields are of particular interest in NMR spectroscopy and 
frequency shift (Knight shift) analysis, especially for sam- 
ples with large magnetic susceptibility.^^ Unlike other 
local field sources, demagnetization fields are always as- 
sociated with the sample geometries. Because the de- 
magnetization field calculation is complex, a general cal- 
culation is rarely available in the literature. Instead, a 
common approach is to simplify the calculation by con- 
sidering more manageable sample geometries like rectan- 
gles, cylinders and spheres and using the average of the 
field as represented by the so called "demagnetization 
factor". ^'^^ The weakness here is that the demagnetiza- 
tion fields are usually not uniform throughout the sample 
even when the sample magnetization is spatially uniform. 

Another approach for the demagnetization field cal- 
culations is to use the fictitious "magnetic charges" as 
studied by Mozurkewich et al.^ Barbara, and Wanas 
et al.^ in which samples with cylindrical geometries 
were again used. 

In this paper, a general calculation for the distribution 
of non-uniform demagnetization fields in paramagnetic 
bulk solids is described. This method involves surface 
currents and volume currents originating from the sam- 
ple magnetization when samples of various geometries are 
exposed to an externally applied magnetic field. Samples 
with cone, ellipsoid, paraboloid and hyperboloid geome- 
tries are investigated. Significant differences in their de- 
magnetization fields are found, revealing the importance 
of the sample geometry for the field. Application is made 
to the theoretical modeling of the ^H-NMR spectra of a 




FIG. 1: (color online) (a) Cartesian coordinate system for 
a sample, where {x,y, z) is the field point inside (or outside) 
the sample volume, and {x' , y' ^z') is the source point (such as 
the point on the surface), (b) The orientation of the applied 
magnetic field Bq. 



two-dimensional (2D) magnetic field-induced supercon- 
ductor (FISC) A-(BETS)2FeCl4 (single crystal), provid- 
ing a good fit to the measured spectra. 



II. DEMAGNETIZATION FIELD 
CALCULATION 

The demagnetization field {B) originates from the 
magnetic polarization of particles. The polarization cre- 
ates magnetic dipoles and forms a sample magnetization 
(M) that leads to a macroscopic surface current density 
(j) on the sample surface and a volume current density 
(J) (if M is not uniform) inside the sample volume. 

Any current produces its own magnetic field in space, 
not excepting for the surface current (density j) and vol- 
ume current (density J) arising from the sample magne- 



tization M, certainly. In any case, it always depends on 
the sample geometry. The demagnetization field B can 
be expressed by 

where /io is the permeability of free space. A' is the sam- 
ple surface area, V is the volume space enclosed by the 
area A', and r is the separation vector, as shown in Fig. 
1 (a), from the source point (x^, , z^) to the field point 
(x, ^, z) in space, i.e., 

f = {x — x')i -\- {y — y')j + (2: — z')k. (2) 

The first term in Eq. (1) is the field produced by the 
surface current density j; it is a point-dependent closed- 
surface integral in a vector space. The second term in Eq. 
(1) is the field produced by the volume current density 
J; it is a point-dependent volume integral in the vector 
space. Here j and J are 

j = Mxh, and J = V x M, (3) 

respectively, where h is the unit vector normal to the 
sample surface, and V is the gradient operator. 

For paramagnetic materials (with non-interacting or 
weakly-interacting moments), the magnetization essen- 
tially has a very accurate linear dependence^^ on the ex- 
ternally applied magnetic field (^o), i-e., 

M = xHo, and ^0 = ^o/mo, (4) 

where Hq is the applied magnetic field intensity, x is the 
material's magnetic susceptibility (isotropic), and Bq = 
Boxi + Boyj + Bo^k^ as shown in Fig. 1 (b). Note, the 
error for the linear relation between M and Hq in Eq. 
(4) is negligible as it is in the order of ^ X^^o, where 
X < ~ 10~^ (cm^/mol-ion) for most known materials. 
Moreover, if magnetic anisotropy is considered, then Mi 
= XijHoj^ where Xij is the tensor element of x {Xij — 
dMi/dHoj) and z, j = x, y, z. 

As seen from Eq. (4), when a sample is placed in a 
homogeneous magnetic field ^o, which is often the case 



as the spatial homogeneity of ^0 is very necessary for 
NMR experiments, there is essentially no spatial vari- 
ation in the magnetization throughout the sample vol- 
ume V . Thus the volume current J vanishes, i.e., the 
2nd term in Eq. (1) J = V x M = 0. In other words, 
the current density only appears on the sample surface 
macroscopically. 

For convenience of practical calculations using Eq. (1), 
we introduce a surface characteristic function f(x'^ y' iz') 
to describe the sample surface where fix'^y'^z') = is 
satisfied [for example, for a unit sphere f{x'^ y' ^ z') = x''^ 

y'^ ^ z'^ — \ \. Then the unit normal h for the surface 
and the surface element dA are 

h = and dA' = j^jy-^dx'dy\ (5) 

respectively, where (Vf)z' is the z'-component of the gra- 
dient oi f{x',y',z'), i.e., (V/)^/ = df/dz'. Note, dA' in 
Eq. (5) comes from dA' = dx'dy' /\n • k\. If n • /c = 0, 
then dA' = dy'dz' /\n • i\ = \{vf) ,\ dy'dz' (if n • z 7^ 0), or 
dA' = dz'dx'/\h • j\ = j^^J^^dz'dx' (if h-j^ 0). 
Thus by combining Eqs. (l)-(5), we have 

= ^ (^"0 X ^) X ^d.'dy'. (6) 

Equation (6) states that if the sample surface is describ- 
able by a surface characteristic function fix'^y'^z')^ then 
the demagnetization field B{x^ y^ z) can be calculated 
straightforwardly with a given external applied magnetic 
field ^0. 

In the following, sample geometries of cone, ellipsoid, 
paraboloid and hyperboloid are examined. 

For example, for a right elliptical cone < z' < c) the 
surface characteristic functions are f{x'^y'^z') = x'^ j €? 
+ y'^j}? - z'^^jc^ and f[x', y' , z') = c - z' , for the lateral 
surface and the top radial surface (at z' = c), respectively. 
Here a, b and c are the semimajor axes along the x'{x), 
y'{y) and z'{z) axes, respectively. With Eq. (6), it gives 
the X, and z components of the demagnetization field 
from the lateral surface current contributions as 



47r 



+ ^) {z - z') - - ^) {y - y') 



x'y' 



TD ^ 



[{x - x')^ + {y- y' f + {z- z' ff^ 
'(-^ - {z - z') - - (x - x') 



x'y' 



B, 



47r jj^.y, 



[{x-x'f + {y-y'Y + {z-z'Yf'^ 

(-^ - ^) {y - y') - + ^) {x - x') 



[{x - x'Y + (y - y'Y + (z - z'Yfl'^ 



—dx'dy' 



-r— -r dx'dy' . 
\z'\ 



(7) 
(8) 
(9) 



2 



Similarly, the contributions from the current on the radial surface [f{x',y',z') = c — z'] are 
I 



Bx = ^ II ^^r^c^a; dy , (10) 

By = ill - — ... — ^^-w, dD 

B, = -^ II """^-^ ' ""^^^ ^,, dx'dy'. (12) 
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Then the total demagnetization field B{x^ z) is the vec- 
tor summation of the two parts, from Eqs. (7)- (9) and 
from Eqs. (10)-(12). 

Similar expressions to Eqs. (7) - (12) can be obtained 
for the ellipsoids, paraboloids and hyperboloids, or any 
other sample geometries as well by using their corre- 
sponding surface characteristic functions fix'^y'^z'). 

Figure 2 shows the calculated result for the demagne- 
tization field magnitudes and directions for samples with 
hyperboloid and ellipsoid geometries with Bq along the 
2;-axis (see Fig. 1). For simplicity the calculations were 
made by setting a = b = c here. The data plots on the 
left side show the magnitude only (plotted as H/Mz vs 
x/a, where is the sample magnetization along the 
z-axis). The vector plots on the right side include both 
the magnitudes and the directions. The corresponding 
results for samples with paraboloid and cone geometries 
are shown in Fig. 3. 

For comparison, the corresponding conditions in the 
calculation for ^o, the a, h and c parameters for all the 
sample geometries, and the scales for the corresponding 
axes of all the plots in Figs. 2-3 are kept the same. We 
also define c/a as the sample aspect ratio. 

Figures 2-3 indicate that the demagnetization fields 
are non-uniform both inside and outside the samples. 
They also indicate that both the magnitude and the dis- 
tribution of the demagnetization fields significantly de- 
pend on the sample geometries at a given Bq (along z). 
The data show that the demagnetization fields are gen- 
erally smaller for the points closer to the z-axis, which is 
understandable even though the cone geometry is slightly 
different in the region of z/a < 0.5 (note, not all the data 
points shown here are close enough to the apex point or 
the lateral surface). The fields are also smaller along 
the direction of Bq since there is no current density (j 
= M X n = 0) at the radial surface (when Bq is aligned 
along z). There is a sharp discontinuity at the edges {x/a 
= 1) except for the hyperboloid, and the integration di- 
verges at the apex for the cone. Among them the smallest 
distribution (a general spatial variation) of the demagne- 
tization fields inside the samples is in the hyperboloid 
under the similar parameter space. 




Demagnetization Fieid Vector 

(Hyperboloid) 



J-t i.t J-t i-t t-L i^-i^.^^^ 



0.2 OA 0.6 0.8 i.O 1.2 lA 1.6 1.8 2.0 

3 r 1 1 \ ^ \ \ ■ 




Demagnetization Field Vector 
(Ellipsoid} 

UJULLU J ^ , 

4 f f n . 

- 1 \ \ f fU:^^„. 



FIG. 2: (color online) Calculated demagnetization field mag- 
nitudes (data plots on the left side) and directions (vector 
plots on the right side) for samples with hyperboloid and ellip- 
soid geometries (shown in the middle), with applied magnetic 
field Bq along the z-axis. 



For further investigation of the field distribution 
throughout the samples, we applied Eq. (6) for indi- 
viduals with fixed geometries. 

As an example. Fig. 4 shows the calculated result for 
the distribution of the demagnetization field Bz{x^y^z) 
[the z-component of B{x^ y, z) (along the direction of 
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FIG. 3: (color online) Calculated demagnetization field mag- 
nitudes (data plots on the left side) and directions (vector 
plots on the right side) for samples with paraboloid and cone 
geometries (shown in the middle) , with applied magnetic field 
Bo along the 2;- axis. 



^o)] inside the samples with various sample aspect ratios 
c/a under a cone geometry. The z-component Bz{x^ z) 
is of interest because it is the component that is along 
the external magnetic field direction that actually con- 
tributes to the NMR Hamiltonian (a negative dot prod- 
uct of the magnetic moment and external field) ii^ii^ The 
line shapes are calculated from Eqs. (7)- (9) and use a 
Gaussian of width (FWHM) 2V2ln2 S, where S = VTO 
(arb. unit). The vertical scales are normalized for com- 
parison. Similar results can be obtained using Eq. (6) 
for samples with ellipsoid, paraboloid, and hyperboloid 
geometries (not shown here). 

Interestingly, Fig. 4 shows that in viewing along the 
applied magnetic field Bq direction, the demagnetization 
field distribution (corresponding to the NMR absorption 
spectrum) is the narrowest and has essentially no differ- 
ence for samples either very thin (small c/a) or very thick 
(large c/a; long sample). This observation is also true 
for samples with a cylindrical geometry, suggesting its 
suitability for samples with other geometries as well. But 
it does have a strong dependence on the sample aspect 
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FIG. 4: (color online) Calculated demagnetization field distri- 
bution for samples with a cone geometry for various sample 
aspect ratios c/a, with applied magnetic field Bo along the 
z-axis {Bo / / c). 



ratio c/a. 

When the sample aspect ratio is small {c/a < 1), the 
demagnetization field variation in space along the direc- 
tion of the applied magnetic field Bq is narrow, while the 
major field variation actually happens most strongly in 
the direction that is perpendicular to Bq. This can be 
seen from Figs. (2)- (3) (case c/a = 1), where the demag- 
netization field variation along x (the radial direction) is 
apparently larger (see the red curves) than that along z 
(view along z from z/a = Otoz/a = l). This is also true 
for all other geometries. 

As the sample aspect ratio c/a increases. Fig. 4 shows 
that the demagnetization field distribution inside the 
sample gets broader with the appearance of a left shoul- 
der (i.e, the line is broadened). This happens starting at 
c/a ~ 1. The shoulder is significant in the range c/a ^ 
3 — 12. As c/a further increases, the shoulder is dimin- 
ished gradually while it gets broader and broader, and it 
becomes insignificant at c/a > ~ 20. However, there is 
little change in the major peak, and the line essentially 
becomes a single narrow peak again. 



III. DISCUSSION 

In this section, we have more comparison discussions 
regarding the sample geometries used in the calculations, 
and we also discuss the result of the application for the 
theoretical modeling of an NMR spectrum. 

Figure 5 shows a typical comparison for the calculated 
demagnetization field magnitude for samples with cone, 
ellipsoid, paraboloid and hyperboloid geometries in their 
middle planes at z/a = 0.5 {c/a = 1), with applied mag- 
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FIG. 5: (color online) Calculated demagnetization field mag- 
nitudes for samples with cone, ellipsoid, paraboloid and hy- 
perboloid geometries in the middle plane z/a = 0.5, with ap- 
plied magnetic field Bo along the z-axis. Here the sample 
aspect ratio c/a = 1. 



FIG. 6: (color online) Calculated distributions of local mag- 
netic fields at a typical proton (HI) sitei^ in a single crystal 
of A-(BETS)2FeCl4 at 20 K and Bo = 9 T. The dashed (red) 
line is for the demagnezation field, the solid (green) line is for 
the Fe^^ dipole field, and the solid (blue) curve shows their 
normalized total. 



netic field Bq along the z-axis. The data comes from a set 
of data shown in Figs. 2-3 (the green curves). Figure 5 
indicates that among them the demagnetization field in 
the cone is the largest, and the smallest demagnetization 
field is in the hyperboloid. Samples with a hyperboloid 
geometry also have the smallest distribution of demagne- 
tization field under similar sample aspect ratio c/a. 

In comparison with other geometries, such as the cylin- 
drical ones,^^'^^ with similar sample aspect ratio {c/a = 
1, for example) and the same Bq [including the field align- 
ment] the demagnetization field magnitude inside a cylin- 
drical sample is also smaller than that for samples with 
a cone geometry correspondingly (like that in a z/a = 
0.5 plane). But the demagnetization field distribution in 
cylindrical geometries is apparently larger and it is also 
more sensitive to the sample aspect ratio c/a}^ A more 
detailed comparison can be made directly by the calcu- 
lation using the similar forms as Eq. (6), i.e.. 
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For samples with a cubic geometry, we made similar 
calculations using Eqs. (13) - (14) and found that their 
demagnetization field magnitudes at the z/a = 0.5 plane 
were ~ 10% larger than those of the cone (here a is the 
cube side length) with the same Bq^ while the field dis- 
tribution inside the cube was the narrowest among all. 

Figure 6 shows the calculated distribution of the de- 
magnezation field using Eqs. (13)- (14) as compared with 
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FIG. 7: (color online) Comparison of the calculated ^H-NMR 
absorption spectrum with the measured spectrum of a single 
crystal of A-(BETS)2FeCl4 at 20 K, for magnetic field Bo = 
9 T aligned with = +85° and = +5°. The dashed (red) 
line is the resull4^ that did not include the distribution of 
demagnetization field contribution at the proton sites, while 
the solid (blue) line is the result that does it. 



the distribution of Fe^+ dipole field at a typical proton 
site (HI site)^in a single crystal of A-(BETS)2FeCl4, 
considering the materials magnetic susceptibility and the 
needle-like sample used for the NMR experiment (at 20 
K and at magnetic field Bq = 9 T). The sample size is 
(1.2 ± 0.1) X (0.065 ± 0.010) X (0.018 ± 0.005) mm^. 
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The material A-(BETS)2FeCl4 is paramagnetic (at 20 
K and 9 T) and has a large magnetic susceptibility x {x 
~ 0.15 emu/mol.Oe) due to the high spin state of the 3d 
Fe^+ ion electron moments (spin S = 5/2).^^ The local 
magnetic fields and the field distributions are dominant 
by the dipolar field from the Fe^+ ions, but they also 
have the contributions from the demagnetization field, 
Lorentz field, exchange field, proton-proton dipolar field, 
hyperfine couplings fields, etc.. Not all of these local 
field contributions are negligible. There are 16 different 
proton sites in each unit cell in this material. In fact, at 
some proton sites the magnitudes of the demagnetization 
fields are even larger than the dipolar fields of the Fe^+ 
ions. This is also true for the Lorentz fields. Thus, the 
demagnetization field in this material could have a signif- 
icant contribution to the local magnetic field magnitude 
and maybe distributions as well, depending on the sam- 
ple geometry and the Bq alignment. 

Figure 7 shows the new result for the theoretical mod- 
eling of a NMR spectrum as compared with the measured 
spectrum of the single crystal of A-(BETS)2FeCl4 at 20 
K, for magnetic field Bo= 9 T aligned with 6 = +85° and 
(j) = +5°, where and (j) describe the direction of Bq rela- 
tive to the crystal lattice. Earlier report did not include 
the contribution of the distribution of the demagnization 
field for simplicity, while the new result does which has 
an application using Eq. (6) [surface currents appearing 
on all the surfaces of the sample] similar to that shown 
in Fig. 6 for all the protons sites in the crystal lattice. 

Noticeably, even though the result shown in Fig. 7 
indicates that the effect of the distribution of demagneti- 
zation field to the shape of the NMR spectrum is not sig- 
nificant in this case, but there is still significant ^H-NMR 
resonance frequency shift contributed by the demagneti- 
zation field. The shift is - 45 x lO-'^Tx 42.5759 MHz/T 
^ 191.6 kHz (i.e., - 0.5%) at - 50 K [note, it is a sig- 
nificant shift even if an NMR frequency shift (Larmor) 
is 0.1%]^^ The small differences between the calculated 
and measured spectra could mainly come from the ir-d 
interactions and d-d interactions of the electron moments 
in the material, which are not included in Eq. (6). 

On an atomic scale, there could be large spatial local 
field variations, thus Eq. (6) can be regarded as the av- 



erage field over a macroscopic small volume dV^ which 
is large enough to contain many atoms (molecules) or a 
unit cell. Thus, in the spectrum modeling, the difference 
of demagnetization field at each proton site in a unit cell 
of the crystal lattice is not considered. 

We have used paramagnetic samples for the demagne- 
tization field calculations, but Eq. (6) can also be used 
for diamagnetic and ferromagnetic materials, as long as 
the linear isotropic relation M = x^o in Eq. (4) ap- 
plies (non-interacting or weakly interacting moments). 
For diamagnetic materials, x < 0? while for ferromag- 
netic and paramagnetic materials X > 0- 
IV. CONCLUSIONS 

In this study, a general calculation of the non-uniform 
demagnetization fields in paramagnetic bulk solids is de- 
scribed and the fields for various sample geometries are 
calculated. The calculations show significant differences 
in their non-uniform demagnetization fields, revealing the 
importance of the sample geometry to the field. Viewing 
from the externally applied magnetic field Bq direction 
when samples are either very thin (sample aspect ratio 
c/a < ~ 0.1) or very long {c/a > ^ 20), we observe that 
the demagnetization field distribution is narrow and in 
this case its contribution to the NMR absorption spectra 
line shape could be neglected. However, the demagneti- 
zation field magnitude and its contribution to the NMR 
frequency shift remain significant, depending on the sam- 
ple magnetic susceptibility as evidenced by our applica- 
tion for the theoretical modeling of the NMR spectra of a 
single crystal of 2D magnetic field-induced superconduc- 
tor A-(BETS)2FeCl4. We suggest that this observation 
be used for all geometries. 
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